Superpenetration of a high energy QQ bound state 
through random color fields 
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The transmission amplitude of a color dipole through a random external color field is computed in 
the eikonal approximation in order to study the absorption of high energy quarkonium by nuclear 
target. It is shown that the internal color state of the dipole becomes randomized and all possible 
color states are eventually equi-partitioned, while the probability of finding a color singlet bound 
state attenuates not exponentially, but inversely proportional to the distance L of the random field 
zone which the dipole has traveled. 



PACS numbers: 11.80.La, 13.85.-t, 24.85. +p 

1. The suppression of the quarkonium production in 
nucleus-nucleus {AB) collisions has been extensively 
studied both experimentally and theoretically since 
it has been proposed as a signal of the QCD plasma for- 
mation. Q The pattern of the observed J/V' suppression 
in light ion induced reactions {AB with B <^ A) ||] as 
well as in pA collisions Q have been well reproduced by 
the nuclear absorption model. [^,|| In this model the ra- 
tio of the observed production cross section ctab in AB 
collisions to the one scaled from the elementary pp colli- 
sions tJpp is given by the simple formula, 
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where L is the effective length of the nuclear medium the 
formative J /i}) (or cc) travels through, (Tabs 6-7 mb 
is the absorption cross section of the " J/?/;" due to the 
collision with individual nucleon, and = 0.16 fm~'^ is 
the mean nucleon density in nucleus. Further suppres- 
sion observed in the Pb-Pb experiments at CERN away 
from this " base-line" has been investigated as an anomaly 
possibly related to the plasma formation. 

The above exponential form of the i-dependence is 
obtained if one assumes a classical stochastic process 
for the multiple collisions of cc pair in nuclear matter. 
One should however note that this treatment is valid 
only when the characteristic time of individual collision 
is shorter than the time interval between the successive 
collisions. It is known that the coherence in the multiple 
collisions changes the high energy asymptotic behavior 
of the scattering qualitatively. | ]To| , pd] | 

Each collision of a cc dipole with nucleon have a char- 
acteristic time scale over which the quantum coherence 
of the wavefunction is important: in the rest frame of 
the pair it may be given by Tc ~ 1/ ls.E ~ 0.3 fm with 
binding energy Ai?, while for moving pairs one should 
take into account the Lorentz time dilatation effect. In 
the Pb-Pb experiment at CERN-SPS, the Lorentz time 
dilatation factor of the J /ijj produced in the midrapid- 
ity region is 7 ~ 10 in the nuclear rest frame, hence 
the cc pair propagates over = ^ 3 fm typically 



in the target rest frame which is of the same order of 
magnitude with the mean internucleon distance in nuclei 
(d = 2 (3/47rno)i/3 ^ 2 fm). In the RHIC (LHC) ex- 
periments, for which 7 ~ 100 (2000) this coherence time 
becomes 30 (600) fm which is considerably larger than 
the mean nucleon distance and become even greater than 
the size of the nucleus. This problem has been pointed 
out recently also in [p^ . 

This simple estimate suggests that the independent 
stochastic collision picture is only marginally satisfied at 
SPS energy and should break down at higher collider en- 
ergies where the coherence effect would become essential; 
it is more natural to consider that the multiple interac- 
tions with nucleons along its path take place simultane- 
ously. It is the purpose of the present work to investi- 
gate how L-dependence of the nuclear absorption is mod- 
ified at collider energies from the naive form (1) due to 
the quantum coherence. For simplicity, here we will not 
go into the discussion of the production mechanism of 
the QQ pair in the nuclear collision which involves color 
octet, as well as color singlet, initial states as produced 
by fusion or fragmentation of incident partons (gluons) in 
nuclear target [ p^p^|jT^ ; we will instead consider a much 
simpler problem how a high energy (color singlet) QQ 
bound state injected in nuclear target will be absorbed 
by multiple collision in nucleus. 

2. Similar problem has been studied previously for the 
absorption of high energy positronium (e"'"e~) while pass- 
ing through metallic foil. ||l^;|l6| It has been found that 
when the positronium energy becomes ultrarelativistic 
E 3> 2r7ie the positronium flux attenuates not by the 
exponential law but by the power law ^ L^^ . This phe- 
nomenon was called superpenetration. Jlq] 

It is easy to see how this happens. [ |l6[ The transmis- 
sion amplitude of the positronium through metallic foil 
may be computed by the eikonal approximation assum- 
ing that the individual members of the positronium travel 
along a straight line trajectory in random atomic electric 
fields. The eikonal phase factor for a particle with charge 
e passing at impact parameter b under the influence of 



the atomic electric field described by the electric poten- 
tial </>(x) is given by ejv J dzip{h, z) hence the transmis- 
sion amplitude of the high energy electron-positron pair 
{v ~ 1), passing at bi = b -I- d/2 and at b2 = b — d/2 is 
simply 



with the small L expansion behavior of the exponential 
decay form (1), if we identify c(r^) — 2a. 

However, for large value of L, P^g{L) deviates from 
the exponential form and exhibits power law behavior, 
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where d = bi — b2 is the transverse size of the dipole and 
q = e J c?zV_L0(b,z) = — e J c?zEj_(b, z) is the net rela- 
tive transverse momentum acquired by the dipole under 
the action of the transverse electric field Ej^ . In this case 
the net momentum transfered to the dipole vanishes due 
to the charge neutrality of the dipole. 

Since the configuration of the electric field seen by 
the pair is random, the relative transfered momentum 
Q = (91 J 92) should be taken as a random variable which 
is to be averaged over with a normalized distribution 
/(q; L); viewing the process as diffusion in the transverse 
momentum space, we expect a simple two-dimensional 
gaussian distribution 



f{q;L)dq 



1 
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with the variance increasing linearly with L: (g^) = cuaL 
where ua is the atomic density and c is a dimensionless 
constant. The penetration probability of the positronium 

state |(^o) is thus given by 



P^„{L) = \{^o\U{q)\ipoW 



dq/(q;i)|F(q)|2 , 



where 



F(q) = J dr_Ldz|^o(r_L,z)pe 



2p-iqr-L 



(4) 



(5) 



is the transverse form factor of the positronium state. In- 
serting (^) into (^) and then performing the integral over 
the transfered momentum q, we obtain an alternative 
expression for P^p„ (L) : 



(6) 

where po{v±) = J dz\ipo{r±, z)\'^ is the probability distri- 
bution of the dipole size r_L in the bound state (i5o(r_L, z). 

For a thin target with small (q^), we can expand the 
exponential in (^) and find 

P<PoiL) ^ 1 - l{rl){q') = 1 - lirDcuAL . (7) 

This result can be obtained by the first Born approxima- 
tion to the eikonal transmission amplitude and coincides 



(8) 



where A = (l/27r) / dq\F{q)\^ = 27r / drj_pl{rj_). Thus 
the probability that the positronium penetrates through 
the metallic foil of thickness L attenuates inversely pro- 
portional to L, not in the exponential form as in (1). We 
observe from the formula that this L^^ dependence 
originates from the depletion of the central value of the 
distribution /(O; L) due to the diffusion of transfered mo- 
mentum in the phase space as L increase. The formula 
(^) indicates also that this is caused by the destructive 
interference of the transmission amplitudes for different 
dipole sizes with |rj^ — r^| > l/y/J/pj. 

Similar asymptotic behavior in multiple collision has 
been noted also by several authors in analyses of high 
energy hadronic process in nuclear target . |17H19|| In par- 
ticular we note the work of Hiifner et al. ||l8| who treated 
the problem of the nuclear absorption of a quarkonium 
in an abelian model of random nuclear color fields with a 
very similar analysis as presented above for the positron- 
ium case and found the 1/L asymptotic behavior for the 
bound state survival probability. However, their conclu- 
sion was plagued by the result of pO| in which the au- 
thors claim the stochastic behavior reappears when one 
includes the non-abelian color degrees of freedom. 

The 1/L-dependence was rediscovered recently by 
Baym et al. p9| by consideration of the Lorentz time- 
dilatation effect for the quantum fluctuation of the trans- 
verse size of the projectile together with color neutral- 
ity argument similar to the treatment in p^ ] 
for positronium absorption. In their semi-classical treat- 
ment, our formula (^) is replaced in effect by 



■(T{r±)nAL 



(9) 



which also gives the 1/L asymptotic behavior if the ef- 
fective cross section of the dipole of size rj^ goes as 
(j{r±) cx r^^.Q Although the quantum coherence between 
successive collisions was seemingly ignored in these treat- 
ments, their result suggests that the L^^ attenuation law 



*We note here that in this semi-classical treatment 1/L arises 
because only small dipole components of the wavefunction can 
penetrate through the medium due to their small absorption 
cross section a{r±). This is not the case in our full quantum 
treatment where the dipoles of large size also contributes to 
the penetration probability. 



is rather universal not specific to the abelian nature of 
the positronium problem. In what follows, we extend the 
works of to include non-abelian nature of the color 
interaction and show that the survival probability of the 
cc bound state will indeed attenuate by the law, 
rather than exponential, at large L. 

3. The new features which arises in the case of non- 
Abelian color interaction is associated with the color de- 
grees of freedom in the internal wavefunction of the pair. 
In the SU(3), there are eight different color orientations 
for the octet cc states in addition to the color singlet con- 
figuration. Action of the external color field will cause 
transitions between these 9 different color states.^ 

To study these qualitatively new aspects, we consider 
here the SU(2) version of the color interaction. The exter- 
nal gauge fields are then given by A^(x) = ^ A'^{x)t'^ /2 
where t'^ (a = 1,2,3) are the (2 x 2) Pauli matrices 
and the internal color states of the dipole are labeled 
by the isospin quantum number: for simplicity, we de- 
notes the singlet state by |0) = |0, 0), and three iso-triplet 
states by |2), |3)) ^ (^(|1, 1) -1)), 1) - 

1)), |1,0)); we use Roman letters a,b, - ■ ■ for label- 
ing the three triplet states and Greek letters, a,(3, - ■ ■ for 
all four states including the singlet state. 

The eikonal transmission amplitude of the color dipole 
QQ may be written as 

UQQ{h,d) ^Ve^^ f dx-^^i i^i-^'^)+ia I dx+A^ {h2,x+) 

where we have used the light-cone coordinates^] = 
{t ± z)/V2 with A± = {A° ± yl^)/\/2, the path-ordered 
product along the integration paths on the light-cone 
— 0) is implied by V and 



a;(xi) = y dx+gA--i^^,x+){T^ - t^)I^ , (11) 
(?(xi) = y" da;+.gA«-(xi,a;+)(Tf +T2")/2 . (12) 

Here t\ (t-i) acts on internal isospin (color) space of Q (ff) 
respectively and = {uji,uj2) and 9 are (4 x 4) matrices 
which operate on the internal isospin (color) space of the 
dipole. 

Some qualitative insights about how non- Abelian color 
interaction works may be obtained by expanding this am- 



plitude by the Born series: 

[/QQ(b, d) 1 + iu}{h) ■ d + i0{h) 



(13) 



where we have shown only the first Born term explicitly. 
The color matrix u) has non-vanishing matrix elements 
only between the singlet state |0) and the triplet states 
\a) : 

(a|C/QQ(b,d)|0)=*(a|u;(b)|0).d + ... 

^4.Vb I dx+gA''-{h,x+) + --- (14) 

while the matrix 6 gives non- vanishing expectation values 
only between the color triplet states: 

{a\UnQ{h,d)\b)=Sab + i{a\e{h)\b) + --- 



Sab + if- 



abcj dx+gA" {h,x+) + 



(15) 



where the antisymmetric tensor eabc is the structure con- 
stant of the SU(2) group. 

Imagine that a color singlet dipole with frozen trans- 
verse size d is injected into a random color field. It will 
then be transformed to color triplet state by the action of 
the field and therefore the probability that it remains in 
color singlet state, Pq.o, will attenuate with the distance 
L as it travels through the random field. In the first Born 
approximation, it may be given by 



Po,o(d;L)^l-5]|(a|a;(b)|0).d|2 



(16) 



where the average is taken over the random field, or 
equivalently, we may consider the matrix a; as a ran- 
dom variable. Since ViA''(x)Vj.4''(x') = (5ab(5y7(x - x') 
by symmetry, we may write 



(a|wj|0)(0|wj|fe) = 5ij5ablQL 



(17) 



where 70 is a constant of the dimension of [volume] ^ . |^ 
Inserting this into ([l^), we obtain 



Po,o(d;L)~l-37oid2 



(18) 



This result coincides with the leading term of the ex- 
pansion of the naive exponential form ([^) if we interpret 



^Action of tlie external field on the color octet cc state may 
also lead to gluon emission from the dipole. We ignore this 
gluon radiation in the following discussion and leave this effect 
for future study. 

*Here we work on the target rest frame where the pair is 
propagating along the trajectory X- = Q. 



^The value of 70 is a measure of the transverse correlation of 
the color electric fields in the target rest frame and we expect 
70 ~ g^^QCD on dimensional ground. In partonic picture, it 
may be related to the saturation scale Qs of the transverse 
parton momentum distribution in a boosted frame. Here 
we take this constant as a pure phenomenological parameter 
which may be determined by comparing our result (|^ to 
data. 



nucleon 




target nucleus 

FIG. 1. Schematic picture of color dipole propagating in 
nucleus which is viewed as a source of random field. 



Cd — 37od^/no as the absorption cross section of the 
dipole of size d. The proportionality of the absorption 
cross section to the square of the dipole size is the man- 
ifestation of the color transparency. pl| , p2| 

If the injected color singlet state is described by the 
internal wavefunction \ipo), we should replace by its 
expectation value (d^)^^ : this procedure may be justified 
by multiplying (a|a;(b)|0) by / dddzip*^{d, z)ipo{d, z) and 
use the closure sum J2 V'n(d, z)Lpn{(i' , z') = S{d—d')S{z— 
z') over all intermediate states n. It then follows that 
if we interpret a^„ = 37o(d^)^„/no as the absorption 
cross section of the bound state, our result coincides with 
the exponential form (|l|) which is based on the classical 
stochastic assumption for the multiple scattering process. 
We will come back to this result later. 

This is a rather surprising result since in our derivation 
we have taken the coherence of scatterings from different 
part of the field fully into account. It is the random field 
average in which eliminates the interference terms of 
these summed amplitude. The formula (|l^), however, is 
based on the first Born approximation and is valid only 
when anoL <C 1. 

4. For large L, we divide the volume occupied by the 
random field into n different small zones along the tra- 
jectory of the dipole in the same spirit as in pO| , p^ . Each 
zone may be considered as belonging to each individual 
hadron (nucleon) in the nucleus pO| | or one could even 
divide the color field in the same hadron into smaller 
portions of cells further ^J]. What is assumed here is 
that the fields belonging to the different zones are un- 
correlated and therefore one can take average over the 
field in each zone independently. The size of the zone I 
characterizes the coherence of the fluctuating color field 
in the nucleus or in the hadron. The random color field 
has been also introduced as a model of the Weizacker- 
Williams fields for small x parton distribution in [p5| . 

The eikonal amplitude is now factorized into the prod- 
uct of the transmission amplitude through each zone: 

UQQ{d,h) = Un---U2Ui , (19) 

where the transmission amplitude for the ith zone, Ui, is 



given by the formula similar to (|1C 

If. ^ gii..><''(b)-d-|-ie<')(b) 

Let us introduce the matrix 



(20) 



■A^i:l';/3,^'(d,d') = (a|C/.(d)|/3)(/3'|f//(d')|a') , (21) 

where the average over the random field in the ith zone 
is implied by the overline. Then the survival probability 
of the color singlet state \(po) after penetrating through 
n zones of uncorrelated random color fields is given by 



= I dddd'po{d)po{d')Ko,o{d,d';n) , (22) 

where po{d) — J dz\ipo{d, z)p is the probability distribu- 
tion of the color singlet dipole of size d in the wavefunc- 
tion ipo{d, z)^^ and the kernel of integration is expressed 
in term of A4 as 

ifo.o(d,d» 



= (0|C/QQ(d,7i)|0)(0|[/J^(d',n)|0) 

where the sum is taken over all intermediate color states 
(ai, = 0, 1, 2, 3) . The formula (^ is non-abelian ex- 
tension of the formula (^ for the positronium penetra- 
tion probability. It still remains to compute the kernel 
function Kofi{d,d';n). 

For a sufficiently small zone the matrix Mi may be 
evaluated by the perturbation theory by expanding 

Ui{h,d) = l + iu}^'\h) ■d + iO'^'^h) 
-i (u;«(b) •d + 6'W(b) 

Calculation of the matrix elements of Mi requires the 
color average of the matrix elements of the products of 
and 6''^*' in addition to (^7|) (with L replaced by I). 
Since there is no preferred color direction and fields in 
the different zones are uncorrelated, we have 



(a|wW|6)(6'|6'«|a'> = 



(24) 



while 



{a\e(^)\b){b'\0U)\a') = eabcea'b'cS^MW , (25) 



** Here we take non-relativistic limit of the light-cone wave- 
function and ignore the mixing of the higher Fock space com- 
ponents which would give energy dependence of the dipole 
distribution. 



where 77 is a constant with the dimension of [length] 
which is also related to the transverse correlation length 
of the gauge fields as 70 and r] ^ g^AgcD, I is the thick- 
ness of each zone. We also need 



(a a; 



(26) 



to compute all matrix elements of the 16 x 16 matrix 
M{d,d') to the order g'^. 

Due to the superselection rules encoded in ([l7|), (p4|), 
( p5| ) and (|2^), the matrix 7W(d,d') becomes block diag- 
onal when it is expressed in the bases of the eigenstates 
of the total isospin and its z components. Hereafter we 
adopt this new representation where Greek letters a, /?, 
stand for (T,r,) = (0, 0), (1, 1), (1, 0), (1, -1). In 
particular, the block matrix for a — a' and /3 — f3': 



(d,d') 



(27) 



is given explicitly by 



M{d, d') 



with 



1 - 3r 
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r' 


r' 
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1 — r — s 
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1 - r - 2s 
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1 — r — s 



(28) 



i7oZ(d^ + d'"), r' = 70M • d' and s = 77(b)/ 
where each column from top to bottom and each row 
from left to right now corresponds to the isospin states, 
(0,0), (1,1), (1,0), (1,-1), respectively. All matrix ele- 
ments with a = a' and (3 ^ (3' which couple the block 
matrix Ma.pid, d') to the other components of the matrix 
A^(*^ vanish by the superselection rules and SapSa'/}' = 0. 
We can express K^xiid^ d'; n) in terms of the the eigenval- 
ues \i and the corresponding (normalized) eigenvectors 
v^'^ of the matrix A/(d, d') as: 



ifo,o(d,d';n) ^Y^^^^v^^f = E ^"/o;' 



(29) 



where /o;i 



,(i)\2 



is the fraction of the color singlet 



component in the ith eigenstate. We find from the ex- 
pression ( ^8| ) of M (d, d') only two eigenstates have non- 
zero mixing of color singlet components for which 



Ai = 1 



+ - ^dX + dt- d\cPL j 70? , 



X2 = l- -[d\ + d 



dX + dt 



d\dl_ 



lol (30) 



and 



, I dl+dl 

/0;1 — 7: ^ 



, I dl+dl 

/0;2 — 7: + 



AjdX + dt~d\d'l 



(31) 



where d| = (d± d')^. Inserting <^ and (jS^) into (|2|), 
we obtain 



i^o,o(d,d';n) 



A? + A5 (4 + rf2_)(A5'-A5) 



AjdX + dt-d\d'^_ 



(32) 



The diagonal components of the kernel A''o,o(d, d; n) 
can be interpreted by the construction as the probabil- 
ity that a color singlet rigid dipolc of size d is found as 
color singlet after going through n uncorrclatcd random 
color zones. Similarly, at d = d' the matrix A/(d, d') re- 
duces to a matrix P(d) = Af(d, d) whose matrix element 
Pa,l3{d) — \ {a\Ui\f3)\^ , can be interpreted as the probabil- 
ity of the dipole of the initial color state a being found in 
the color state (3 in the final state after going through the 
ith random field zone. Since P is a symmetric stochastic 
matrix P {pap) with J^aPap = Y^pPaP = li it l^^s 
an eigenvalue of 1 for the eigenvector v = (1,1,1,1)^ 
and other eigenvalues are smaller than 1. Moreover, the 
matrix P" will approach with increasing n to the ma- 
trix with all elements equal to 1/4. This latter property 
implies that the transition probability from the initial 
isospin (color) singlet state to all color states are equal; 
namely, there is equi-partitioning of all different isospin 
states in the final state. The physical reason for the ap- 
pearance of the stochastic color equi-partitioning is the 
cancellation of all interference terms due to the random 
color averaging. 

The approach to the color equilibrium can be calcu- 
lated explicitly from (32) as 



ifo,o(d, d; n) = + \{l 4r)" = \ + ^e^'^ (33) 

where k = — ln(l — 4r(d, d)) ~ 4r = 470/d^. 

This stochastic behavior of the motion in the internal 
color space was already noticed in It resembles 

to the exponential decay if one ignores the asymptotic 
value 1/4 (1/9 in the case of SU(3)). This is not the 
entire story, however: the relaxation process in the color 
space is limited by the finite number of color degrees of 
freedom and the bound state survival probability still 
attenuates by the 1/L law due to the other continuous 
degrees of freedom, namely the transverse size d of the 
color dipole which is frozen during the collision, as we 
shall show now. 

The survival probability of the bound state can be ob- 
tained by inserting (^2|) into ( p2|) and performing the in- 
tegral over d and d'. It is important here to take into 
account the interferences of the transmission amplitudes 
at different dipole size d 7^ d' as we have seen in the cal- 
culation of the positronium penetration probability: the 
information about the internal momentum transfer q by 
the scattering in the random field is contained in the de- 
pendence on d — d' of the kernel Ko,o{d,d';n), namely 
in its off-diagonal components. 



We first check that our formula contains the thin target 
result by setting n = I: 



= 1 - 3(d2)7oL = 1 - L/L 



l--{dl + dl)jol 



(34) 



where we have used L = I and Labs = l/3(d^)7o may be 
interpreted as the absorption length of the bound state. 
This is exactly what we have obtained already. 

For thick targets with large n, the terms with the 
largest eigenvalue, Ai, dominate so that 



rfdrfd'po(d)po(d') 



AJdi+dt-dldl 



A? . (35) 



To estimate the remaining integrals we observe that the 



integrand is symmetric by the interchange ^4 
equivalently by d' <-> —d', and the factor 



d-, or 



A? 



n „n In Ai 



- ^ {d'^+-\-d'i--^dl^dt~d^^d^_)'yonl 



(36) 



becomes 1 when d = ±d' but decreases very rapidly as 
d ± d' increases. Expanding the exponent in terms of 
d?. = (d — d')^ and setting d = d' in the prefactor, we 
may estimate roughly as 

S^M ^ 2 I dddd'pUd)\e-i(^-<''^'-<"^ ^ ^ , (37) 



where the constant Lq is given by 

Lo = 2xix2x A(d^)Labs 
4 



(38) 



with A = 27r J (idp§(d) and we have used 3(d^)7o = auQ. 
This behavior is essentially the same as we have seen in 
the positronium case, (||), except for the factors 2 x i x 2: 
the first factor 2 arises from the symmetry by d+ <-> 
d- which is specific to the case of SU(2) color charge; 
the second factor j originates from the equi-partitioning 
in the color degrees of freedom and will be replaced by 
= 1/9 for SU(3) color; finally the last factor 2 has 
arisen because only one of the two eigenvalues Ai and 
A2 contributes at large L. In the case of SU(2) color 
charge these factors cancel accidentally, while in the case 
of SU(3) color charge, they are replaced bylx-p-x2 = |. 

Having established the asymptotic behavior of the sur- 
vival probability S^pg (L) we examine how the transition 
from the thin target result ( p^ ) to the thick target asymp- 
totic behavior ( p7| ) takes place using the gaussian normal- 
ized dipole size distribution: po(d) = -Kje~'^ /'^o which 
may be obtained from the ground state wave function of 



the non-rclativistic heavy quark of the reduced particle 
mass /X in the harmonic oscillator potential with the fre- 
quency Lo = l/pdQ. In this case, A(d^) = 1 and two 
length scales iabs and Lq coincide. The survival proba- 
bility S^pg then becomes a function of a single dimension- 
less variable x = L/Li^i^g- The numerical result is plotted 
in Fig. ^ and compared with the exponential decay form 
obtained by the exponentiation of the thin target result 
(H): 1 - L/Labs e'^l^'^^'^. It is seen that the expo- 
nential absorption formula gives an overestimate of 20 % 
for the suppression at L = Labs in the case of SU(2) color 
charge, while this is slightly reduced in the case of SU(3) 
color charge as shown in Fig. 3. 

In summary, we have studied the absorption of a high 
energy QQ bound state passing through random color 
fields taking fully into account the quantum coherence 
of the multiple scatterings while assuming the transverse 
size of the dipole is frozen by the Lorentz time dilata- 
tion. It was shown that the absorption is weaker than 
that given by the exponential damping form commonly 
used in phenomenological models for nuclear absorption 
and is given instead by the power law inversely propor- 
tional to the thickness of the target asymptotically. Al- 
though our calculation of penetration probability of a QQ 
bound state in nuclear target is not directly applicable to 
the quarkonium production problem in nuclear collision, 
our result suggests that a special caution is needed to use 
the naive nuclear absorption model of quarkonium sup- 
pression at high energies, especially at RHIC and LHC 
energies. However, it remains to be seen how the coher- 
ence effect will show up in the nuclear collisions when 
one includes the production mechanism of QQ pair [Q. 
We are now working on the problem and the result will 
be reported elsewhere. 
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